Exploring the graphene edges with coherent electron focusing 
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We study theoretically the coherent electron focusing in graphene nanoribbons. Using semiclas- 
sical and numerical tight binding calculations we show that armchair edges give rise to equidistant 
peaks in the focusing spectrum. In the case of zigzag edges at low magnetic fields one can also observe 
focusing peaks but with increasing magnetic field a more complex interference structure emerges 
in the spectrum. This difference in the spectra can be observed even if the zigzag edge undergoes 
structural reconstruction. Therefore transverse electron focusing can help in the identification and 
characterization of the edge structure of graphene samples. 
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I. INTRODUCTION 

Transverse electron focusing (TEF) is a versatile ex- 
perimental technique which has been used in metals to 
study the shape of the Fermi surface and the scatter- 
ing on various surfaces and interfaces^. The accessibility 
of the quantum ballistic transport regime in two dimen- 
sional electron gas (2DEG) opened up the way to the 
experimental demonstration of the coherent electron fo- 
cusing 2 in GaAs heterostructures as well. 

The geometry of the coherent electron focusing Ref. [H 
is shown in FigJT] The current is injected into the sam- 
ple at a quantum point contact called injector (I) in per- 
pendicular magnetic field. If the magnetic field is an 
integer multiple of a focusing field Bf ocus , electrons in- 
jected within a small angle around the perpendicular di- 
rection to the sample edge can be focused onto a sec- 
ond quantum point contact (the collector, denoted by C 
in Fig. [IJ which acts as a voltage probe. Therefore, if 
the collector voltage is plotted as a function of magnetic 
field one can observe equidistant peaks at magnetic fields 
B = p* Bf ocus (p — 1, 2, 3, ... ). The first focusing peak 
corresponds to electrons reaching the collector directly, 
i.e. without bouncing off the edge (see Fig. [Tj. Subse- 
quent peaks correspond to trajectories bouncing off the 
edge p — 1 times before reaching the collector and there- 
fore their presence can attest to the specular nature of the 
scattering at the edge. In other words, the p > 2 focus- 
ing peaks can give information on the scattering process 
taking place at the edge of the sample. 

The edge structure of graphene 3 nanoribbons 4 ^— and 
graphene flakes^ have recently attracted a lot of inter- 
est because it strongly influences the nanoribbons' and 
flakes' electronic and magnetic properties ^ 13 i 15 . Theo- 
retically, the most often studied edge structures are the 
armchair and zigzag ones which have recently been ob- 
served experimentally as welli&£. Density functional cal- 
culations suggested that other types of edges might also 
be present, comprising pentagons and heptagons of car- 
bon atoms- 13 . Experimental evidence for this type of edge 



reconstruction has indeed been found very recently 14 . 
The effect of the hydrogen concentration of the environ- 
ment on the edges has also been studied and further pos- 
sible edge structures identified^. Experimentally how- 
ever the identification and characterisation of the edge 
structure has often been a challenge 6 . 




FIG. f : (Color online) Schematic geometry of the transverse 
electron focusing setup. A graphene nanoribbon is contacted 
by an injector (I) and a collector (C) probe and perpendicular 
(to the graphene sheet) magnetic field is applied. Classical 
quasiparticle trajectories leaving from the injector at normal 
direction, depending on the magnetic field, can be focused 
onto the collector. 

In this paper we show that coherent electron focusing 
can be used in ballistic graphene samples to study the 
properties of the edge structure. We argue that in the 
case of armchair edges one would see equidistant peaks in 
the focusing spectrum at integer multiples of a focusing 
field Bf ocus . In contrast, for zigzag and reconstructed 
zigzag 13 (reczag) edges only the first few focusing peak 
would be identifiable and for stronger magnetic fields a 
more complex interference structure would appear in the 
focusing spectrum. The presence or absence of focusing 
peaks at stronger magnetic fields can therefore discrimi- 
nate between armchair and zigzag (reczag) edges. 

Our main results are summarized in Fig. [5J Us- 
ing the tight-binding model for graphene with nearest- 
neighbour hopping we numerically calculated the trans- 
mission probability T(B) from the injector to the collec- 
tor as a function of the magnetic field B for three dif- 
ferent types of graphene nanoribbons: armchair, zigzag 
and zi gzag with reconstructed edges, denoted by zz(57) 
in Ref.[l3l As Fig. [5] shows one can indeed observe peaks 
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FIG. 2: The results of tight binding calculations for the trans- 
mission probability T(B) as a function of magnetic field, (a) 
for an armchair; (b) for a zigzag nanoribbon. The black verti- 
cal lines indicate the positions of the focusing peaks predicted 
by the semiclassical theory, see Sections IPT1 and Hill 



and zigzag nanoribbons making use of the semiclassical 
theory of graphene, introduced in Refs. [20l|2ll 

The rest of the paper is organized as follows. First, 
Section [TT] we discuss the boundary conditions for the 
semiclassical theory on the edges of armchair and zigzag 
graphene nanoribbons and derive the quantization con- 
dition for edge states in magnetic field. The quantization 
condition then allows us to calculate the band structure. 
Using this in Section ITTT1 we show that the focusing spec- 
tra of armchair and zigzag nanoribbons are different in 
strong magnetic fields. We then turn to the comparison 
of the results of numerical calculations and the theoreti- 
cal predictions. In Section Hvl we discuss some details of 
the numerical calculations which underpin the theoreti- 
cal approach presented in Sections [TT] and Mil Finally, in 
Section [V] we give our conclusions. 



II. SEMICLASSICAL THEORY OF EDGE 
STATES IN GRAPHENE 



in the focusing spectrum of these nanoribbons. (In the 
one orbital per site approximation we used in the com- 
putations the focusing spectrum of the reczag nanorib- 
bon is very similar to the simple zigzag one's therefore 
we only show the latter here.) In the case of the re- 
constructed edge we took into account that the hopping 
between the atoms is different on the heptagons and pen- 
tagons than in the bulk of the graphene. To obtain re- 
alistic nearest-neighbour hoppings we employed ab-initio 
calculations. (See Section [TV] for the details of the ab- 
initio method and Ref. [l6| for the actual parameters of 
the calculations). We consider the case where there is a 
finite carrier density in the sample, i.e. the Fermi energy 
is well above of the Dirac point (this ensures that the 
semiclassical approach we take in Sections [TTJ and IIIII is 
justified). We will focus on the transmission peaks which 
can be observed for B/Bf ocus > 1 because they can give 
local information on one of the edges (see Section [TTJ) . 
They should also be present if e.g. a graphene flake is 
contacted by two probes, as long as the edge between 
the probes is not disordered. The first peak in the trans- 
mission for all three types of nanoribbons can be found 
at B/B focus « 1 where B focus = 2M£. (k F is the Fermi 
wavenumber (measured from the K point) and L is the 
distance between the injector and the collector). While 
for armchair edge there are well-defined peaks whenever 
B is integer multiple of Bf ocus [see Fig. Ufa)], in the case 
of zigzag edges (both the ideal and the reconstructed 
one) a more complex interference pattern emerges for 
B/Bf ocus > 6 showing many oscillations but not a clear 
peak structure [Fig. HJb)]. These results imply that a) 
one can distinguish the armchair and zigzag edges by 
their high magnetic field focusing spectrum; b) the dif- 
ference in the focusing spectra persists even if the zigzag 
edge undergoes structural reconstruction. We will ex- 
plain the differences in the focusing spectra of armchair 



We start our discussion by establishing the link be- 
tween the theory of semiclassical approximations for 
graphene and the boundary conditions for Dirac fermions 
on honeycomb lattice. If the magnetic length Is = 
\J H/\eB\ is much larger than the lattice constant of 
graphene, the general energy-independent boundary con- 
dition has the form of a local restriction on the compo- 
nents of the wave function ^ at the edge (£) 18 ' 19 . It 
can be cast into the following form: M = W where the 
4x4 matrix M may be chosen as Hermitian and unitary 
matrix: M — M\ and M 2 = I. One can show-i 9 - that 
demanding: a) that the probability current normal to the 
boundary be zero; b) that the boundary should preserve 
the electron-hole symmetry of the bulk; c) and finally, 
assuming that the boundary conditions do not break the 
time reversal symmetry, leads to the following form of 
the matrix M: M = (vt ® ncr), where a — (cr x , cr y , a z ), 
t = (jx ,T y ,T z ) and <7j , Tj are Pauli matrices acting in the 
sublattice and valley space, respectively. Furthermore, u 
and n are three dimensional unit vectors, restricted to 
two classes: zigzag-like (y = ±z, n = z, where z is the 
unit vector perpendicular to the plane of the graphene 
sheet) and armchair-like (y z = n z = 0). Additionally in 
both classes n _!_ n^, where ri£ is a unit vector in the 
plane of the graphene sheet and it is perpendicular to 
the edge. 

Let us now consider the reflection from an edge of a 
nanoribbon in more details. It follows from the form of 
the boundary conditions described earlier that the wave- 
functions \tg on the boundary is proportional to the 
eigenvectors Z^ corresponding to the doubly degener- 
ate unit eigenvalue of the matrix M. In other words, 
*f = n ± Z ± e lkx where *+ (#") and Z + (Z~) corre- 
spond to isospin vector u (-v), 77 ^ are amplitudes and 
k x is the wavevector component along the (translation- 
ally invariant) edge. On the other hand, one can show 
that in magnetic fields where Ib 3> Af = 2n/kp (Xp 
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is the Fermi wavelength) the wavefunctions can also 
be written as a superposition of an incident \ n and re- 
flected ^out plane wave: \& f = ^f n + r ± ^E , Q Ut where 
are reflection amplitudes. By equating the two forms of 
H!^ the coefficients f ± and r? ± can be easily obtained for 
any boundary condition described by the matrix M. For 
instance, in the case of armchair edge, where the isospin 
vector can be parametrized as v = (cost/?, sin </>,0) and 
n = (1,0, 0) T , the eigenvectors of M with unit eigenval- 
ues are 



a) 



The ansatz for tyf can also be written as 



m 



1 



1 



(2) 



where a is the incidence angle (measured from the x axis, 
see Fig. [T]). Equating the two forms of i.e. VE^ = 
i] ± Z ± e lkx and Eq. ([2|) a straightforward calculation gives 



= Exp(iA$^); where A<E>+ = a, A$ a = a + it. 

(3) 

In the case of zigzag edges the boundary can be char- 
acterized by a superlattice vector T = ma\ + la^ where 
m ^ I are integers and a±, a,2 are the lattice vec- 
tors of graphene 19 . The boundary condition is given 
by M — £ m i (t z (g) o z ) where £ m j = sgn(m - I), where 
sgn(. . . ) is the sign function. After performing analogous 
calculations as for the armchair edge one finds that 

ff = Exp(iA$J); where A$J = a - tt =f £ mi a<. (4) 

The physical meaning of Eqs. ©, ([3]) and (Q} is the 
following: assuming specular reflection in classical pic- 
ture the momentum of the particle is rotated upon re- 
flection. Since the quasiparticles in graphene are chiral, 
the change in the direction of their momentum rotates 
the pseudospin as well and therefore it leads to a change 
in the phase of their wavefunction. As one can see from 
Eqs. ([3]) and (0| this phase shift depends on the type 
(armchair v. zigzag) of the edge as well. This will be 
important when we use semiclassical quantization to ob- 
tain the band structure of graphene armchair and zigzag 
nanoribbons. 

A second ingredient in the calculation of the edge 
states is the semiclassical theory introduced in Refs. [2CJ, 
SH. As it has been shown^i Refs. HJU one 
can introduce a classical Hamiltonian H e (p,r) = 

v f\J (p — eA) 2 + V(r) for graphene to describe the clas- 
sical motion of electron-like quasiparticles. Here p = 
(pxiPy) is the canonical momentum, A(r) is the vector 



potential describing any external magnetic field and V(r) 
is scalar potential which is taken to be zero throughout 
our discussion, V(r) = 0. If the edges of the graphene 
ribbon do not break the translational invariance the sys- 
tem is integrable because the longitudinal momentum p x 
and the energy are conserved (see Fig. [T] for the choice 
of the coordinate system). Using the Landau gauge 
A(y) = (By, 0,0) T for the vector potential to describe 
perpendicular magnetic field pointing into the —z direc- 
tion, the Hamiltonian T-L e will not depend on the x co- 
ordinate and the only nontrivial quantization condition 
is related to the motion perpendicular to the graphene 
edges (in the y direction). In general it can be written as 
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S y + 7 + A$f = 2tt 



(5) 



where S y = £ p y dy is the classical action and the inte- 
gration is over one period of the motion perpendicular 
to the edges. Furthermore, A$£ are the phase shifts 
coming from the reflections at the edges; is a pos- 
itive integer, fi is the Maslov index counting the num- 
ber of caustic points and finally 7 is a Berry-phase-like 
quantit y 20 ' 21 which can be calculated from the equation 
^7 (r(i)) = g (V r x 4tt) . Here the classical trajectory 



r(t) of quasiparticles is given b y " = V p H e (p, r). 

One can show that in perpendicular magnetic field 7(1") 
is half of the deflection angle of the momentum of the 
quasiparticles. 

In Landau gauge the classical action S y can be calcu- 
lated in the same way as in Refl2ll. The detailed pre- 
sentation of the semiclassical quantization of armchair 
and zigzag nanoribbons in perpendicular magnetic field, 
discussing all the possible classical orbits and the corre- 
sponding 7 and A$g phases is left for a forthcoming pub- 
lication. For our purposes it will be sufficient to consider 
only those classical orbits which correspond to skipping 
motion along the boundaries, see e.g. in FigfTJ If the 
magnetic field is strong enough such that the diameter 
of the cyclotron orbit is smaller than the width of the 
nanoribbon (denoted by W in FigO}, i.e. for 2R C < W 

where R c = is the cyclotron radius at Fermi energy 

Ep, one can show that only these skipping orbits corre- 
spond to current-carrying states. The phase shift A<E>£ 
coming from the reflection at the edge is given by Eq. ^ 
and Eq. Q for armchair and zigzag edges, respectively. 
Considering first armchair nanoribbons, using Eqs. (J3j 
and ([3]) the semiclassical quantization of skipping orbits 
reads 



k F R c 



+ /3± + -sin2/3± 



2w{n ± ± 



(6) 



where kp = 7-7^, fin 1S the angle with the y axis under 
which the cyclotron orbit is reflected from the bound- 
ary, and +(— ) corresponds to isospin u(—v). Moreover, 
n+ = 0-..n+ ax {n~ = l...n" a J where n^ ax is the 
largest integer smaller than ifcj?i? c =p j. This quanti- 
zation conditions holds for both metallic and insulating 
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armchair nanoribbons. The physical meaning of Eq. ([6]) 
is that the enclosed flux by the cyclotron orbit and the 
edge of the nanoribbon is quantized and it must equal 
/i/e(n ± ± 1/4). It is interesting to note that for 2DEG 
with hard wall boundary conditions the quantization con- 
dition for edge states looks very similar, except that one 
would have to omit the ± sign and on the right hand side 
one would have (n + |). 

The semiclassical quantization for edge states in the 
case of zigzag nanoribbon can be obtained in the same 
way as for an armchair nanoribbon, the only difference is 
that instead of Eq. one has to use Eq. (Q} . One finds 
that it is given by 



(7) 



k F R c + - sin 2/?± J -vUl (| + 



= 2ir I - - 
4 



Here ± corresponds to the product f£ m / = ±1, where 
both v and £ m z can take on values ±1. The value of 
v depends on the isospin: it is 1 (—1) if the isospin is 
v = z (y = —z) and £ m / = ±1 is introduced before 
Eq. (j4]). The range of the quantum number n* 1 is the 
same as for the armchair case. One can notice that as 
compared to Eq. ([6]) there is an extra term on the left 
hand side of Eq. . The origin of the difference in the 
focusing spectra of armchair and zigzag nanoribbons, to 
be discussed in Section IIII1 can be traced back to this 
term in the dispersion relation. 

We have calculated the band structure for armchair 
and zigzag nanoribbons in homogeneous perpendicular 
magnetic field both semiclassically, using Eqs. ©, (|7J), 
and numerically in tight binding approximation. As one 



low energies E/fnUc < 1/2. (The dispersionless state of a 
zigzag nanoribbon at zero energy cannot be described by 
semiclassics either.) However, the dispersionless sections 
of the band structure corresponding to Landau levels at 
finite energies can be calculated semiclassicall y 20 ' 21 . (We 
note that similar results have been obtained in Ref. 
using the Dirac-like Hamiltonian for graphene.) 



III. MAGNETIC FOCUSING IN GRAPHENE 
NANORIBBONS 

Having obtained the quantization condition for edge 
states in Eqs. © and |J7|), the calculation of the mag- 
netic fields B focus i where the transmission between the 
injector and collector is sharply peaked, follows the rea- 
soning of Ref. 0- The ballistic transport along the edges 
of a nanoribbon can be understood in terms of the edge 
states described in Section [TTI because they are the prop- 
agating modes of this problem. If the injector is narrow 
(~ \f) one can assume that it excites these modes co- 
herently. Therefore, as long as the distance between the 
injector and the collector is smaller than the phase co- 
herence length, the interference of the edge states can be 
important. For a given Ep the interference of the edge 
states, labelled by n ± in Eqs. © and ([7]), is determined 
by the phase factors exp(ifc^L). Here the wave numbers 
&;„ are given by k~ = /cFsin/3^ and L is the distance 
between the (very narrow) injector and collector. Both 
in the armchair and in the zigzag case one can show that 



around n 



c /2, which corresponds to 



< 1, 



in good approximation the angles /?„ depend linearly on 
ri ± . 

Expanding Eq. © around (3^ — one finds that 
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FIG. 3: Comparison of the results for the band structure ob- 
tained from semiclassical quantization given by Eqs. ((6}, 
(circles) and from numerical tight binding calculations (solid 
lines). The energy is in units of fju c — \/2^-- The magnetic 
field is given by W/Ib = 8.97. (a) armchair nanoribbon; (b) 
zigzag nanoribbon, in the vicinity of the K point. 

can see in Fig.|3]the agreement between the semiclassical 
and the tight-binding calculations is very good except for 



fc„ L 



kpL 



2n d 



(8) 

where C a = jkpL. This result means that if L/2R C is 
an integer (or equivalently, B/Bf ocus is integer, where 



B 



focus 



2hkr. 



r ) some edge channels, with quantum num- 
bers n centred around n^j aa ./2, can constructively in- 
terfere at the collector. Other edge states, to which the 
linear expansion shown in Eq. ((5J) cannot be applied, give 
rise to an additional interference pattern which does not 
have a simple periodicity. The expression for the focus- 
ing field B focus is formally the same as for 2DEG 2 - but 
the dependence of k F on the electron density is different 
in the two systems. 

Repeating the expansion for zigzag nanoribbons using 
Eq. (J7J we find that to linear order 



k„ L 



1 ttL 



it 



Rr 



-c, 



(9) 



where k — Ef/FujJc is the filling factor and C z = 
— jkpLj — 2-r- . For nanoribbons where the cyclotron ra- 

dius becomes smaller than the width of the ribbon such 
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that k 1 is satisfied and therefore 1/ (l =p j^t) w 1 the 
focusing field Bf ocus can be defined in the same way as 
for armchair ribbons and the transmission peaks are at 
the same focusing fields as in the armchair case. How- 
ever, according to Eq. (|9]) as k gets smaller for increasing 
magnetic field (assuming fixed Fermi energy) the focus- 
ing field for the two isospin index becomes slightly dif- 
ferent, moreover, the subsequent focusing fields for each 
isospin index are no longer equidistant. (Examination of 
the higher order terms not shown in Eq. ^ reveals that 
also the number of interfering edge channels becomes dif- 
ferent for the two isospin indices.) As a result, a more 
complex interference pattern is expected to appear than 
in the case of armchair nanoribbon and a simple Bf ocus 
can no longer be defined. 

Turning to the comparison of the numerical results 
shown in Fig. [5] with the analytical predictions given in 
Eqs. ([8} and (|9]), one should focus on the B/Bf ocus > 1 
regime, where the edge states are the current carry- 
ing modes. One can see that the peak positions for 
the armchair nanoribbon are in very good agreement 
with the prediction of the semiclassical theory. Up to 
B/ Bf ocus S 5 the focusing peaks are also clearly dis- 
cernible in the focusing spectra of zigzag nanoribbons 
but for stronger magnetic fields a more complex inter- 
ference pattern emerges. The peak at B/Bf ocus = 7 for 
example can clearly be seen for the armchair case. In con- 
trast, for the zigzag edge a number of oscillations with 
similar amplitudes can only be observed. The difference 
between focusing spectra of the armchair on one hand 
and of the zigzag (and reczag) edges on the other hand 
is even more noticeable in stronger magnetic fields, i.e. 
for B I Bf ocus > 7. (Eventually, at very strong magnetic 
fields the cyclotron radius of the quasiparticles would be- 
come comparable to the width of the collector and the 
system would therefore be equivalent to a Hall bar. We 
will not consider this regime here.) Note, that for a typ- 
ical electron density of n e — 2.5 x lO 16 ^- and assuming 
L = 1/im the focusing field is Bf OC us ~ 0.37T, which is 
experimentally feasible. 

In the case of the calculations presented in Fig. [2] the 
width W of the nanoribbon was larger than L/2 (L is 
the injector-collector distance). This meant that for the 
Fermi energy we used already the p = 1 focusing peak 
could be observed. For narrow ribbons however it can 
easily happen that for given Ep and W the first few 
peaks would not appear, until for sufficiently strong field 
the cyclotron radius R c becomes smaller than W and 
the condition 2 * p * R c = L is met for some integer 
p. For narrow zigzag and reczag ribbons therefore there 
might be no observable focusing peaks because for the 
magnetic fields where R c would become small enough, 
the k » 1 condition [see below Eq. (|9])] is no longer 
satisfied. There could be of course oscillations or peaks 
in the transmission even in this case but those would not 
have the simple periodicity that focusing peaks have. 



IV. DETAILS OF THE NUMERICAL 
CALCULATIONS 

Let us now briefly discuss some of the details of our 
numerical calculations for the transmission probability 
T(B) shown in Fig. [2] The injector and the collector were 
modeled by heavily doped graphene and the transmission 
was calculated employing the Green's function technique 
of Ref. [l?]. The graphene nanoribbons are assumed to be 
perfectly ballistic and infinitely long. This means that 
the left and right ends of the nanoribbons act as drains 
which absorb any particles exiting to the left of right. To 
simulate the effect of finite temperatures we used a sim- 
ple energy averaging procedure in the calculation of the 

transmission curves: T(B) = J T(B, E) (-^TzT^) dE 
where fo(E) is the Fermi function. The actual results 
shown in Fig. [2] were calculated at T = IK temperature. 
As it can be expected, higher temperatures tend to smear 
the curves while at lower ones an additional fine structure 
appears. 

For simple armchair and zigzag nanoribbons we as- 
sumed that the hopping parameter between the atoms is 
the same everywhere on the ribbons. In contrast, for the 
zz(57) edge we took into account that the hopping param- 
eter changes on the pentagons and heptagons at the edges 
with respect to its the bulk value. We calculated the 




FIG. 4: The reconstructed zigzag edge. The numbers indicate 
the change of the hopping parameter on particular bonds with 
respect to its bulk value, 7 = — 3.39eU. 

hoppings with density-functional tight-binding (DFTB) 
method, using the hot bit code^ This approach yields 
the hopping matrix elements as two-center integrals for 
pseudo-atomic orbitals that result from a straightfor- 
ward pathway of ab initio calculations, as described in 
Ref. [23|. The density-functional parametrization yields 
a valid description of the covalent bonding in carbon 
nanomaterialsi2£ Fig U shows the relative hoppings from 
DFTB-optimized geometry, note how the appearance of 
a triple bond in the armrest parts affects also the n- 
electron hopping, due to the significant reduction in bond 
length^ In the one orbital per site approximation we 
employed here the band structures of zigzag and reczag 
nanoribbons, apart from the close vicinity of the Dirac 
point, are very similar (see Fig. [5]). This explains why 
they have very similar focusing spectra as well. 
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FIG. 5: Comparison of the band structures of zigzag (solid 
line) and reczag (circles) nanoribbons from tight binding cal- 
culations at the K point. The strength of the magnetic field 
is given by W/Ib = 9.14. 

V. SUMMARY 

In summary, we studied coherent electron focusing in 
graphene nanoribbons using both exact quantum calcu- 
lations and in semiclassical approximation. We found 
that in the case of armchair edges the transmission peaks 
are at integer multiples of B/Bf ocus whereas for zigzag 



edges such a simple rule holds only when the filling factor 
k is much larger than unity. For zigzag nanoribbons in 
stronger magnetic field, when the filling factor is of the 
order of one, a more complex interference pattern can be 
observed in the transmission and the emergence of this 
interference pattern can be understood from semiclassi- 
cal calculations. The presence of focusing peaks at low 
magnetic fields can therefore attest to the quality of the 
edge structure, while measurements at stronger fields can 
discriminate between armchair and zigzag edges. Our 
numerical calculations on zz(57) edges suggest that the 
above conclusion holds even if the zigzag edge is struc- 
turally reconstructed. Although we considered in our an- 
alytical calculations nanoribbons whose both edges had 
perfect armchair or zigzag structure, we expect that our 
findings are more generally valid. Namely, they rely on 
the properties of edges states localized at one of the edges 
therefore imperfections of the other edge should not af- 
fect the focusing peaks. Therefore this technique should 
be applicable to study the edges of graphene flakes as 
well. Finally, an interesting extension of our work would 
be to consider the focusing spectrum of other proposed 15 
edge types and possibly taking into account more than 
one orbital per site in the transport computations. 
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